Higher Charmonia and X,Y,Z states with Screened Potential 
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We incorporate the color-screening effect due to light quark pair creation into the heavy quark- 
antiquark potential, and investigate the effects of screened potential on the spectrum of higher 
charmonium. We calculate the masses, electromagnetic decays, and El transitions of charmonium 
states in the screened potential model, and propose possible assignments for the newly discov- 
ered charmonium or charmonium- like "X,Y,Z" states. We find the masses of higher charmo- 
nia with screened potential are considerably lower than those with unscreened potential. The 
Xc2(2P) mass agrees well with that of the Z(3930), and the mass of ■!/;(4415) is compatible with 
tp{5S) rather than tp{4:S). In particular, the discovered four Y states in the ISR process, i.e., 
y(4008), y(4260), F(4320/4360), F(4660) may be assigned as the ip{3S), ip{4:S), ip{3D), ipiQS) 
states respectively. The X(3940) and X(4160) found in the double charmonium production in e"'"e~ 
annihilation may be assigned as the ric{3S) and Xco(3P) states. Based on the calculated El transition 
widths for Xci(2P) 'yJ/^ and Xci(2P) 77/)(2S') and other results, we argue that the X(3872) 
may be a Xci(2P) dominated charmonium state with some admixture of the O'^D*" component. 
Possible problems encountered in these assignments and comparisons with other interpretations for 
these X, Y, Z states are discussed in detail. We emphasize that more theoretical and experimental 
investigations are urgently needed to clarify these assignments and other interpretations. 

PACS numbers: 12.39.Jh, 13.20.Gd, 14.40.Gx 



I. INTRODUCTION 



In recent years, a number of charmonium-like states, 
the so-called "X,Y,Z" mesons, have been found at B 
factories and other experiments (for reviews see e.g. (T], 
0, Q). These states all lie above the open charm (e.g. 
D{*)D{*)) thresholds, and the observed masses and de- 
cays make the identifications of these new states very 
puzzling. Apart from the conventional charmonium, 
many exotic candidates, such as molecules, tetra-quarks 
and charmonium- hybrids, are suggested. However, de- 
spite many exciting hints, none of the exotic states has 
been firmly established so far, and more theoretical stud- 
ies are needed to explain the existing data and to confront 
new experimental tests. In the domain of conventional 
charmonium, one of the main difficulties to assign the 
new mesons as excited charmonia is that the observed 
masses do not fit the predictions of potential models 
with linear confinement potential (the quenched poten- 
tial). However, the situation can be more complicated 
than the simple potential model calculations by includ- 
ing the coupled channel effects, or the string breaking 
effects. These effects could make the masses of higher 
charmonia lower than the potential model predictions. 
One distinct example is the state Z(3930) observed in 
the two-photon process [3], which is now identified with 
the 2P charmonium Xc2' but its mass is about 40 McV 
lower than the prediction given by the quenched potential 
model [5]. Another example is the X(3940) observed in 
double charmonium production in e^e^ annihilation 6], 
which is likely (though not necessarily) to be the 35' spin- 
singlet charmonium state 77" Q , but this assignment will 



imply that the mass of the 3S pseudoscalar charmonium 
is smaller than the prediction given by the quenched po- 
tential model by about 100 MeV Q. Therefore, it is 
possible that the quenched potential model may overes- 
timate the masses of charmonia in the energy region well 
above the open charm thresholds. 

Although very successful in the prediction of the char- 
monium spectrum below the open-charm threshold, it is 
well-known that the quenched potential model, which in- 
corporates a Coulomb term at short distances and the 
linear confining interaction at large distances [1, Q , will 
not be reliable in the domain beyond the open-charm 
threshold. This is because the linear potential, which 
is expected to be dominant in this mass region, will be 
screened or softened by the vacuum polarization effects 
of dynamical fermions [10]. Unfortunately, this screening 
effect has not been directly detected with the standard 
Wilson loop technique in the unquenched lattice calcu- 
lations. The reason might be that the Wilson loop op- 
erator almost decouples from the physical ground state 
at a large lattice distance that consists of two disjoint 
strings [ll|. In other words, to simulate the screening 
effect out with the Wilson loop technique needs a very 
long lattice time, which might be too far beyond the abil- 
ity of the simulation at present. On the other hand, the 
screening or the string breaking effect has been demon- 
strated indirectly by the investigation of the mixing of a 
static quark-antiquark string with a static light meson- 
antimeson system in ri/ ~ 2 lattice QCD [l2|. This effect 
is also confirmed by the calculations within some holo- 
graphic QCD models [l^ recently. However, since the 
simulations of lattice QCD still have large uncertainties 
and difficulties in handling higher excited states, in order 
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to investigate the screening effects on the charmonium 
spectrum it should be useful to improve the potential 
model itself to incorporate the screening effect, and com- 
pare the model predictions with experimental data. Such 
screened potential models [13, [H, [H, [i3| were proposed 
many years ago in the study of heavy quarkonium and 
heavy flavor mesons, as well as light hadronsfTl|. In re- 
cent years the screened potential models have again been 
used to investigate the heavy quarkonium spectrum and 
leptonic decay widths [l^. In addition to heavy hadrons, 
the spectra of light hadrons have also been investigated 
with the screened potentialfl^, and it is argued that the 
large degeneracy observed in the excited meson spectrum 
by the Crystal Barrel Collaboration in pro ton-antiproton 
annihilation in the range 1.9-2.4 GeV'21 may indicate 
the flattening of the confinement potential due to the 
color screening effects [20j. On this point, it is important 
to further examine experimentally whether or not the 
linear Regge trajectory of the meson spectrum can hold 
for even higher excited mesons. More complete data are 
needed in the future to clarify the issue of the light meson 
spectrum regarding the color-screening effect. Although 
no definite conclusion can be drawn at present, it is cer- 
tainly useful to study the color-screening effects on the 
mass spectra for both heavy and light hadrons, especially 
for the newly discovered higher excited states. 

The effect of vacuum polarization due to dynamical 
quark pair creation may also be compensated by that of 
the hadron loop induced by virtual D meson pairs in the 
so-called coupled-channel model [1, [13, [H, [2J| , for char- 
monium system. In Ref. [1^, we compared the coupled- 
channel model with the screened potential model in char- 
monium spectrum in the mass region below 4 GeV. With 
the same quenched limit, the two models are found to 
have similar global features. It is not surprising since, in 
the quark-meson duality picture, the two models may de- 
scribe roughly the same effects. In practice, calculations 
with the screened potential model are simpler to deal 
with, though detailed predictions for the spectrum can 
be somewhat different from the coupled-channel model 

In this paper, we calculate the mass spectrum of the 
charmonium especially the higher charmonium using a 
non-relativistic Schrodinger equation with the Coulomb 
potential plus a screened linear potential which is nearly 
the same as that in [13, [3] but with slightly different 
parameters (see Sec. II for details). Spin-dependent in- 
teractions are considered perturbatively. With one ad- 
ditional screening parameter /i, the model predicts that 
the masses of higher charmonium are lowered, compared 
with the quenched linear potential, and this mass sup- 
pression tends to be strengthened when the charmonium 
states vary from lower levels to higher ones. For instance, 
we find that the calculated mass of Xc2 is 3937 Mev, fit- 
ting well the experimental value of Z(3930) [3|, and the 
mass of ip{5S) is close to the observed ^'(4415). Con- 
sequently, this mass spectrum will leave more room for 
new assignments for some of the observed X, Y, Z states 
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FIG. 1; Comparison of the unscreened potential V^(t-) (dashed 
line) and the screened potential Vacr (solid line) for r = 0.1-4 
fm with parameters taken from (|10p . The asymptotic limit of 
Vscr is shown by the dotted line. 



in the charmonium family. These possible assignments 
will be suggested and discussed in detail in this paper. 

In the following, we first introduce the screened po- 
tential model in Sec. II, and then study some decay and 
transition processes of charmonia in Sec. III. In Sec. IV 
we will discuss possible assignments for the observed 
charmonium(-like) states. A summary will be given in 
Sec. V. 



II. THE SCREENED POTENTIAL MODEL 

As a minimal model describing the charmonium spec- 
trum we use a non-relativistic potential model with 
screening effect being considered [13, [iBl ■ We use a po- 
tential as 

Vscr{r) = Vvir) + Vs{r), (1) 
Vvir) = (2) 

Vs{r)^X{ ), (3) 

where is the screening factor which makes the long 
range scalar part of Vscr (?") flat when f ^ ^ and still lin- 
early rising when r <^ —, X is the linear potential slope, 
and ac is the coefficient of the Coulomb potential. Fig- 
ure 1. shows the screened potential departure from linear 
in large length with the parameters given in (jlOp . 

For hyperfinc interactions, we only consider the spin- 
dependent interactions which include three parts as fol- 
lows. 

The spin-spin contact hyperfinc interaction is 

Hss^^^Ur)Sc-S,, (4) 
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where Sa-(r) is usually taken to be S{f) in nonrela- 
tivistic potential models, but here we take Sa{r) ~ 
(ct/a/tt)'^ e"*^ as in Ref. Q since it is an artifact of an 
0(u^/c^) expansion of the T-matrix [2^ in a range com- 
parable to 1/TOc. 

The spin-orbit term and the tensor term take the com- 
mon forms: 



III. SOME DECAY PROCESSES 

A. Leptonic decays 

The electronic decay width of the vector meson is given 
by the Van Royen-Weisskopf formula 29] with QCD ra- 
diative corrections taken into account 30]. 
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{-Vy{r)-Vy{r))T. 



(5) 



These spin-dependent interactions are dealt with pertur- 
batively. They are diagonal in a | J, L, S" > basis with the 
matrix elements 
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T^e{nS) = -r72— l^ns(0)| (1 - ), 



(11) 



T^AnD) ^ J|^|i?::^(0)p(l - IJ^), (12) 

(6) 2M^^m^ 3 TT 

where M„5(M„£)) is the mass for nS{nD), Cc = | is the 
c quark charge in unit of electron charge, a is the fine 
structure constant, i?„5(0) is the radial S wave-function 
at the origin, and i?„^(0) is the second derivative of the 
radial D wave-function at the origin. 

(7) Combined with the parameters (fTO|) , we get the results 
listed in Table HIl 



<L-S>= [J{J L{L + 1) - Sis + l)]/2 (8) 



B. Two-photon decays 



and the tensor operator T has nonvanishing diagonal ma- 
trix elements only between L > spin-triplet states, 
which are 



<3 Lj\T\^Lj <; 



6(2L+3) 



6(2L-1) 

For the model parameters, we take 



J = L + 1 

J = L 

J = L-1. 



(9) 



ac = 0.5007, A = 0.21GeV^, 
fi = 0.0979Gey, a = 1.362GeF, 
mc = 1.4045Gey, as = 0.26, 



(10) 



where ac ~ as{rncVc) and as ~ as{2mc) are essentially 
the strong coupling constants but at different scales. The 
former is for large distances and used to determine the 
spectrum while the latter is for short-distances and used 
for QCD radiative corrections in charmonium decays (see 
below in next section). Here fi is the characteristic scale 
for color screening, and l/fi is about 2 fm, implying 
that at distances larger than l//i the static color source 
in the cc system gradually becomes neutralized by the 
produced light quark pair, and string breaking emerges. 
With these values of the parameters for the potential, we 
can calculate the spectrum of the charmonium system. 
The results are shown in Table |T1 For comparison, we 
also list the experimental values :2T] and those predicted 
by the quenched potential model '5] in Table HI 



In the nonrelativistic limit, the two-photon decay 
widths of ^S'o, '^Pq, and can be written as [3l| 
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(13) 



r^-CPo - 77) = ^^"^^-'^"-("^'^ , (14) 
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(15) 



The first-order QCD radiative corrections to the two- 
photon decay rates can be accounted for as [3l| 
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)] ,(17) 
(18) 



sensitive to the details of potential near the origin. So 
we take 

n'So -> 77) — > r(%p 77) pe.p«(^^^ (^9) 

to eliminate this uncertainty. 

In the nonrelativistic limit, we can also replace m^. by 
M/2, where M is the mass of the corresponding charmo- 
nium state. The results are listed in Table Hill 
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El transitions 



IV. DISCUSSIONS 



For the El transitions within the charmonium system, 
we use the formula of Ref . [s^ : 

rEi(n2S+iLj^^,2S'+iL;j,+^) 

= ^Cf,Sss'ela\{f\r\^)\^E'^ (20) 

where is the emitted photon energy. 
The spatial matrix element 

<f\r\i>= / Rf{r)R,{r)r^dr, (21) 
Jo 

involves the initial and final state radial wave functions, 
and the angular matrix element C/^ is 

2 



Cf^ = max(L, L') (2J' + 1) 



L' J'S 



(22) 



J L 1 

Our results are listed in Table HVl The widths calcu- 
lated by the zeroth-order wave functions are marked by 
SNRq and those by the first-order relativistically cor- 
rected wave functions are marked by SNRi. 

For the first-order relativistic corrections to the 
wave functions, we include the spin-dependent part of 
(HI),®,® and the spin-independent part as ^] 

H,r = -i + iv'H'(r) 
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Pi ■ r- 



(23) 



where Pi and P2 are momenta of c and c quarks in 
the rest frame of charmonium, respectively, which sat- 
isfy Pi = — P2 = P, 3 is the unit second-order tensor, 
and {{ }} is the Gromes's notation 

{{A-n-B}} = ^{AB : n+A-m+B-nl+n : AB), (24) 

where 3? is any second-order tensor. 

Note that we do not include the contributions from 
the scalar potential in Hgj since it is still unclear how to 
deal with the spin-independent corrections arising from 
the scalar potential theoretically. 

The results of non-relativistic linear potential model [5] 
marked by NR are listed for comparison. Our results 
for the El transition widths of SNRq are slightly larger 
than those of the NR's mainly due to the relativistic 
phase space factor of e['^'^''/m|'^'^' involved in the NR's 
in Q , where Ej'^'^'' is the final state charmonium total 
energy, and m|'^'^^ is the initial state charmonium mass. 
Both SNRq and NR's results of El transitions are larger 
than experimental values, but we see that in SNRi the 
predicted widths get decreased and fit the experimental 
values rather well as long as the first-order relativistic 
corrections to the wave functions are included. 



A. Z(3930) 

The Z(3930) was discovered by the Belle Collabora- 
tion [J as an enhancement in the DD invariant mass 
near 3.93 GeV/c^ in the 77 collision with the statisti- 
cal significance of 5.3(7. Results for the mass, width, and 
product of the two-photon decay width times the branch- 
ing fraction to DD are: 

M(Z(3930)) = 3929 ± 5 ± 2 MeV, (25) 
r(Z(3930)) = 29 ±10 ±2 MeV, (26) 
T-y-yB{Z {S9S0) ^ DD) = 0.18 ± 0.05 ± 0.03 KeV, (27) 

respectively. The production rate and the angular distri- 
bution in the 77 center-of-mass frame suggest that this 
state is the previously unobserved Xc2 (2P) Q . 

The mass of Xc2(2P) in our model is 3937 MeV, which 
is consistent with the experimental value in (j25J) , whereas 
the mass in the quenched potential model [5| is larger 
than the experimental one by 40-50 MeV (see Tabcl[I]). 
As we have mentioned, this is one of our motivations 
to reexamine the charmonium spectrum in the screened 
potential model. 

The open-charmed decays of Xc2(2P) were studied in 
Ref. P (the mass was set to be 3931 MeV), and the 
total width is predicted to be 35 MeV with the branch- 
ing ratio S(Z(3930) DD) ~ 74%. Together with the 
data in (P?]) . one can get the two-photon decay width 
of Z(3930) as about 0.16-0.33 KeV, which is consistent 
with the predicted value 0.23 KeV for Xc2(2P) in our 
mode (see TabelUlB. 



B. V(3770), i/'(4040), i/)(4160), i/'(4415) 

Before the discovery of X{3872), there are only 
four well-established charmonium states above the DD 
threshold. They are -0(3770), ^(4040), ^(4160) and 
■!/'(4415). They all have quantum number J^^ = 1 . 
Conventionally, they were assigned as mainly '0(1'^ Z3i), 
tpiS^Si), il){2^Di) and V'(4^5'i), respectively. These as- 
signments are consistent with the predictions given by the 
quenched potential models. However, in our screened po- 
tential model, the -0(4415) should be assigned as -0(5'^ S":) 
(see Table HI). 

The experimental di-electron width of 0(3770) is 
0.259 ± 0.016 KeV, larger than that expected for a pure 
D-wave state, which is probably due to mixing with the 
2-^ 5*1 state induced by the coupled channel effects as 
well as by the tensor force. In the following we assume 
-0 = V(3686) and = 0(3770) to be admixture of 
2'^iS'i — l^Di charmonium states and 



10') = |23Si)cos6l+ |l3Di)sin6', 
IV-") = -|23Si)sin6l+ |l^Di>cos6', (29) 



(28) 
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where 9 is the mixing angle. This angle can be estimated 
by comparing the experimental values of the di-electron 
widths of V'(3770) and ■!/'(3686) with theoretical predic- 
tions for pure l^Di and 2^5*1 states (see Table HH). The 
result is given by 



-12° or 6'w25°. 



(30) 



Thereinto, 6 — 25° is not compatible with the width of 
V"' Xco7, so we take 9 = -12°. Note that 9 « -12° is 
consistent with the results of Refs. [H, [H, [H, Ull . 

The S-D mixing may be more serious for i/)(4160) (with 
ree(V'(4160)) = 0.83± 0.07 KeV^), if one assign it 
as the 2^Di state, since its observed di-electron width 
is comparable to that of ^"(4040) (with ree(-0(4O4O)) = 
0.86± 0.07 KeV[23). Extracted from the ratio of di- 
electron widths of ■!/'(4160) and "0(4040), the mixing angle 
between 2^Z)i and 3"^ 5*1 can be as large as —37° in our 
model and —35° in Refs. [i^, H^. Moreover, there could 
also be mixing between 2D — AS states. This indicates 
that the observed di-electron widths for higher charmo- 
nia (above the open charm threshold) can be altered by 
the S-D mixing effects due to coupling to the decay chan- 
nels. Moreover, their masses can also be modified by the 
S-D mixing. 

The open-charmed decays of these states are not stud- 
ied in this paper. However, one may expect that the 
corresponding calculations and results in our model will 
be similar to those in the quenched potential model. In 
Ref. 01 , the authors evaluated the open-charmed decay 
widths of these states using the '^Pq quark pair creation 
model, and the results are consistent with experimental 
measurements both for the total widths and the decay 
patterns. Especially, they predicted that the main de- 
cay modes of ■;/;(4415) (as V'(45') in their model) are DDi 
and DD2 ^, and the DD2 mode was confirmed by Belle's 
measurement [ill recently. The lesson from both the the- 
oretical calculation and the experimental measurement is 
that the higher excited charmonia tend to decay into ex- 
cited charm mesons rather than the S-wave charm meson 
pairs. This might be due to relativistic suppression and 
the node structures of the wave functions of the higher 
excited states. 



C. y(4008),F(4260),y(4320/4360),F(4660) and 
X(4630) 

The vector state F(4260) was first discovered by 
the BaBar Collaboration ,28] as a relative narrow 
peak around 4260 MeV in the J/ip tt+tt" distribution 
in the initial state radiation (ISR) process e+e^ 
7iSR>^/'0 TT+TT". This state was also observed by 
CLEO [11] and Belle [49], and there is also a relative 



broad structure near 4.05 GeV, the so-called F(4008), in 
the Belle data. Recently, BaBar [50|] updated their mea- 
surement and gave the mass and width of F(4260) which 
are consistent with Belle's measurements [49|. However, 
Babar has not yet seen the structure around 4.05 GeV. 
The mass and width of the broad structure Y(4008) are 
given bvjiol] 

M(y(4008)) = 4008 ± 401^^'' MeV, (31) 
r(y(4008)) ^ 226 ± 44 ± 87 MeV. (32) 

The averaged mass and width of Y(4260) are given bv[27j 

M(r(4260)) = 4263t^ MeV, (33) 



r(y(4260)) = 95 ±14 MeV. 



(34) 



BaBar also found a structure around 4.32 GeV [5l| in 
the ISR process e^e~ 7iSR^'7r~''7r~ with 

M(y (4320)) = 4324 ±24 MeV, (35) 
r(y(4320)) = 172 ±33 MeV, (36) 

while in the same process the Belle Collaboration ob- 
served two relative narrow peaks around 4.35 GeV with 

M(r (4360)) = 4361±9±9 MeV, (37) 
r(r(4360)) = 74 ± 15 ± 10 MeV, (38) 



and 4.66 GeV \m with 



M(y(4660)) = 4664 ± 11 ±5 MeV, 
r(y(4660)) = 48 ±15 ±3 MeV. 



(39) 
(40) 



^ In this paper, inclusion of the charge conjugate mode is always 
implied. 



Aside from the broad structure y(4008), which might 
be related to the ip{iS) state ?/'(4040) and even the ip{2D) 
state 7/;(4160), the other three F-states are considered 
to be difficult to assign as conventional charmonia since 
their masses are inconsistent with those predicted by the 
quenched potential models [H^l , and there are even no 
enough unassigned states in the charmonium spectrum 
to accommodate them. As a consequence, these F-states 
are interpreted totally or partly as exotic states, such as 
ccg hybrid [54] 1 cqcq tetra-quark state [S^, [H, Wh : t>ay- 
onium [H^ and molecule [5^. As the only exception, in 
Ref.jeSI the Y(4260) was interpreted as the ipiAS) char- 
monium. Recently, in Ref. 5^ the authors assigned 
i:{3D) and ^(55*) to the ^(4325/4360) and F(4660) 
states, although the predicted masses are higher than 
the experimental values by 50 to 100 MeV. 

The situation changes greatly in our screened poten- 
tial model. In the earlier calculations fldi. [iTj . the mass of 
-0(45') was just around 4260 MeV, the mass of Y(4260). 
Since the predicted higher chamonium spectrum is com- 
pressed in the screened potential model, there is enough 
space to accommodate these y-states. Specifically, in 
the model of this paper the masses of ip{4:S), i^{3D) and 
■0(65') (see Table H]) are predicted to be 4273, 4317 and 
4608 MeV, which are roughly compatible with the ob- 
served masses of y(4260), r(4325/4360) and r(4660). 
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respectively. The small mass discrepancies between the- 
oretical predictions and experimental data may be either 
due to the experimental errors or, more likely, due to the 
theoretical uncertainties, especially the complicated S-D 
mixing effects, such as the mixing among 45", 31? and 55* 
states. 

The di-electron widths of the pure S'^Di and G'^Si 
states in our model are 0.97, 0.044 and 0.49 and KeV, 
respectively. Experience about the large S-D mixing (es- 
pecially between 7/;(4040) and ^'(4160)) discussed in the 
last subsection tells us that the large S-D mixing may 
change these di-electron widths to moderate values. As- 
suming the di-electron widths of r(4260), ^(4325/4360) 
and y(4660) to be all about 0.4 to 0.5 KeV, one can ex- 
tract the branching ratios B{Y —^JPi^tD 7r'''7r~) from 
the experimental measurements [49l |52|. They are 
all about 1-2% and the corresponding widths r{Y 
il;{tp') TT+TT") = 1-2 MeV. These widths look too large 
compared with r(i/'' J/^ tt+tt") ~ 100 KeV 
However, once the higher charmonium lies well above the 
open charm threshold, the di-pion transition rate may 
be enhanced dramatically by final state interactions be- 
tween charmed mesons, which are produced in the decay 
of the charmonium. This is similar to the case in the 
T{5S) di-pion transitions [HH, where the rescattering ef- 
fects between B^*") — B^*") mesons are expected to play a 
crucial role in enhancing the di-pion transition rates of 
T{5S) into T(1S', 25*). Furthermore, if the broad struc- 
ture y(4008) is indeed due to ^'(4040) as well as V'(4160), 
we will have a clear experimental hint for the large di- 
pion transition rates of 3S as well as 2D charmonium, 
which are well above the DD threshold. 

It is useful to emphasize that for all the four newly 
discovered Y states the measurements [H lislligl. ligllsol. 

m 

ree{Y)xB{Y^^7r+n-)^Oil0eV), (41) 

where tp means J/^p or ^(2S), and the measured total 
widths of about 0(100 MeV) altogether may imply that 
they have similar properties. Therefore a coherent inter- 
pretation for these four Y states is needed. A likely expla- 
nation is that they are conventional charmonium states, 
though other interpretations are also possible and even 
more interesting. Note that in the hybrid scenario if one 
of y states, e.g., the Y(4260) is the 1 hybrid, then one 
needs to understand why the others, which can no longer 
be accommodated as 1 hybrids in this mass region, 
have similar properties to the Y(4260). On the other 
hand, there are some considered difficulties to assign Y 
states as conventional charmonia, and the most serious 
one to assign Y(4260) as the 4S-doniinated charmonium 
seems to be the observed dip rather than a peak in the R 
value scanned in e'^e~ annihilation [27J around Y(4260) 
(this difficulty is common to all resonance interpretations 
of Y(4260)). A possible explanation for the dip is the 
destructive interference between the continuum and the 
resonance. If without any resonance in this region, the 
continuum contribution should be generally smooth. An- 



other difficulty is the nonobservation (not a peak but a 
dip observed) of the decay modes DD, DD*, D*D* at 
the Y(4260) (for a recent report by BaBar see Ref . [6^ ) . 
The dip in the observed charmed meson pairs is proba- 
bly related to the dip in the R value, since the latter is 
the measurement of the hadron production cross section 
in e^e~ annihilation. The dip in R is just the reflection 
of the dip in the resonance decays to hadrons (only the 
charmed hadrons are relevant here). They may all be 
caused by the interference effects. Moreover, the above 
mentioned difficulties are not only for the charmonium 
assignment but also for other interpretations of Y(4260). 
One needs to understand the dip in R around Y(4260) 
if one tries to interpret Y(4260) as a resonance no mat- 
ter which kind of resonance it is. Nevertheless, it is in- 
structive to search for decay modes involving the P-wave 
charm mesons e.g. DiD and other higher charm mesons, 
apart from the S-wave charm meson pairs, because the 
higher charmonium may prefer decays to higher charm 
mesons or multi-mesons, due to the form factor suppres- 
sion with higher momentum released in decays to lower 
charm mesons, and also due to the node structure of 
higher charmonium state. In this regard, we note that a 
main decay mode of the ■0(4415) is D2D (47| . 

Very recently, Belle reported a new vector state 
X(4630) [1^ which was found as a threshold enhance- 
ment in the A^A~ distribution in the ISR process 
e'^e^ — > 7isrA+A^. The mass and width are fitted to 
be 

M(X(4630)) = 4634t^+| MeV, (42) 
r(X(4630)) = 92tfit2i MeV, (43) 

which are roughly in agreement with those of y(4660). 
Assuming that AT (4630) is the same state as the y(4660) 
{'ip{6S) in our model) with the di-electron width of about 

0. 5 KeV, one can extract the partial width 

r(A(4630) ^ A+A") - 10 MeV, (44) 

1. e., the branching ratio of about 10%. Such a large bary- 
onic decay width certainly deserves further studying. 

To sum up for the discussion in this subsection, our 
assignments for these newly discovered Y states appear 
to be consistent with the 4S, 3D, 6S charmonium mass 
spectrum predicted by the screened potential model, and 
other properties may also be understood in these char- 
monium interpretations. But the issue is far from being 
conclusive, and many theoretical and experimental inves- 
tigations are apparently needed to clarify these assign- 
ments with other more interesting interpretations such 
as hybrids and tetraquarks. 

D. X(3872) 

The A (3872) was first observed by Belle ^ in 
the J/ip ■K^i:~ invariant mass distribution in B'^ — s- 
J/tp TT'^TT~ decay as a very narrow peak (Tx < 2.3 
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MeV) around 3872 MeV. The mass of X(3872) in the 
J/tp TT+TT" mode was recently updated by CDF Collab- 
oration [HI] as 

M(X(3872)) = 3871.61 ± 0.16 ± 0.19 MeV, (45) 

which is very close to the £)°Z)*° threshold m(D"L)*") = 
3871.81 ± 0.36 MeV jlll. The spectrum of the di-pion 
indicates that they come from the p resonance and the 
charge parity of X is ev en [ 64l| . Moreover, analyzes both 
by Belle [67| and CDF (68| favor the quantum number 
J^^ = 1++. 

The product branching ratio B{B+ XK+) ■ B{X 
J/^p{-K+TT-)) is about 7-10 X 10"^ [Hill. With the 
rate of this mode, the relative rates of other decay modes 
of X(3872) areinllzl 



which implies that X(3872) is an isoscalar. The most 
serious problem of the molecular model, in our opin- 
ion, is that it is difficult for a loosely bound state to 
radiatively transit into exited charmonium, such as ip' ^ 
through quark annihilation or other mechanisms. Model 
calculations [t^] predict the ratio 



^(^ - V/ 7) ^ 4 , 



whereas the experimental value of this ratio [72| is 



i?;^^/^^ = 3.4±1.4. 



(52) 



(53) 



Rip' ^ 



= 1.0±0.5, 



^ 0.33 ±0.12, 



1.1 ± 0.4. 



(46) 
(47) 
(48) 



Most of the above problems for the molecular model 
can be resolved if one can assign X(3872) as a conven- 
tional charmonium. As a J^'-' = 1++ state, the only 
candidate is the x'd whose mass is about 3.90 GeV in our 
model (see Table. I). The 30 MeV difference between the 
predicted mass and the experimental one in (j45p can be 
further reduced if the coupled channel effects are taken 
into account [1^. It is the S-wave coupling of x'd to 
]jO£)*o -j-jjg^-j. i^gjijig to lower the mass of x'd towards the 
threshold of D^D*'^ . This is related to the cusp effect at 
the D°D*° threshold [79^. 

The charmonium candidates of X(3872) were sug- 
gested [sol ? ] soon after it was found. However, 
these suggestions were almost given up, after the isospin- 
violating decay X J/ipp was confirmed. Because of 
the coupled channel effects, the x'd will mix with nearby 
opened D'^D*'^ component. Such a mixed charmonium 
model for X(3872) was proposed in Ref.^1] and Ref.[82]. 
M(X(3875)) = 3872.6t'^ l ± 0.4 MeV, (49) Differing from the molecular models, the D°D*° compo- 



B{X- 


^ J/ip uj) 


B{X- 


. J/i; p) 


B{X- 


. J/^ 7) 


B{X - 


. Jl^ p) 


B(X- 


V' 7) 


B{X- 


> J/^ p) 



It is interesting that another narrow structure was 
found in the D°D°tt° ^ or the D^D*° [^ invari- 
ant mass spectrum near 3875 MeV, which is a lit- 
tle higher than that in in the decays iS+Z^ 
D°D*°{D°n°)K+^". Recently, Bellelz^ updated the 
measurement on X(3875) and improved their fitting 
method and found 

5-1-0.5 



which is consistent with that in ((45|) . Provided that the 
two X-states are the same, one can extract the ratio 



Rdd* 



B{X D^D*^) 
B{X ^ J/ij p) 



9±2 



(50) 



from the Belle data [75| . 

The X(3872) is widely accepted as a molecule can- 
didate of D^D*^ in S-wave [3, [ttI since its mass is 
very close to the D^D*^ threshold. This assignment can 
also give a natural explanation of the J^'~^ of X(3872) 
and predict the ratio R^u) ~ 1 [ttI ]. which is in agree- 
ment with that in (|46p . However, as a loosely bound 
state of D^D*^, it should be difficult to be produced 
in B-decays or in p — p collision at the Tevatron. For 
example, a model calculation [t^I shows that in the 
B'^ decay a molecule X(3872) has a branching ratio 
B{B+ XK+) = (0.07-1.0) X 10-*, whereas the ex- 
perimental rate tends to exceed this upper limit. Fur- 
thermore, Belle Collaboration has observed X(3872) in 
the neutral channel B° -> X{J/ip TT+n-)K° with 5.9(7 
significance and with a rate almost as large as that of the 
charged channel [70l | 



B{B° X{3872) K°) 
B{B+ X{3872) K+) 



0.82 ±0.22 ±0.05, (51) 



nent mixed in the 1"*""*" charmonium is just a hadronic de- 
scription for effects of the vacuum polarization induced 
by the dynamical quark pair creation and annihilation. 
Thus, the mixed chamonium is as compact as the con- 
ventional charmonium. As a result, the production rates 
of X(3872) should be large and equal in both the neutral 
and charged channels in B meson decays The pro- 
duction rate of X(3872) in p—p collisions at the Tevatron 
should also be large, comparable to that of Xci(l^') (but 
somewhat reduced due to a smaller cc norm in the mixed 
charmonium model of X(3872)). 

Using the final state rescattering mechanism, one 
may explain the isospin violating decay X(3872) 
J 14' P |83l |. The isospin violation, which is implied by the 
ratio Rtfii^ in ([46]), is expected to be mainly due to the 
difference between the thresholds of D^D*^ and D~^D*~ , 
and the larger phase space of J/4 P than that of J/tp lu 
also favors the J/tp P decay [83 . ISSj . In addition, the 
ratio Rdd* in ((50|) may also be accounted for provided 
that the X(3872) lies below or just a little amount above 
the D°D*" threshold 

The El transition rates of x'd sensitive to the rela- 
tivistic corrections due to the node in the 2P wave func- 
tion, especially for the one x'd ^ J 14 7- In our model, 
after relativistic corrections are taken into account, the 
transition widths T{x'd ^ J/4{4') l) = 45(60) KeV 
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fsee Table II V|) . The corresponding ratio i?^/^/^ ^ ~ 1.33, 
which is much larger than the one predicted by the molec- 
ular model in ([52|) and in rough agreement with the ex- 
perimental value (|53p . Different treatments or different 
parameters in the relativistic corrections can result in 
very different estimations for the rate of x'ci ~* J /''P 7 
(2P-1S transition), while the rate of x'ci ^ "0' 7 (2P- 
2S transition) can only be changed a little. For exam- 
ple, Ref. [13 gives T{x'ci ^ J/^{-<P') l) = 11(64) KeV, 
and the corresponding ratio R^ij/^ ^ ~ 6. Thus, in 
the mixed charmonium model for X(3872), the expected 
range of the ratio may be 



R^ij/^^ = 1.3-6.0. (54) 



If we use the calculated T{x'ci i^'l) = 60 KeV as in- 
put for X(3872) ^'7, and use the experimental results 
(|46ll . (|47| . (|48ll . and ([5^1), as weU as the width of decay to 
light hadrons (assuming T{xci{'^P) ^ light hadrons) ~ 
r(xci(lf ) light hadrons) « 600 KeV), we wih get 
the total width of X(3872) to be about 1400 ± 300 KeV, 
which is compatible with the measurement (it can be fur- 
ther reduced when the cc norm in X(3872) is less than 
one). 

The nature of X(3872) can also be uncovered by the 
pole structure of the scattering amplitude involving the 
resonance near the D^D*'^ threshold. This study is also 
needed to explain the different peak locations of X(3872) 
in the JMtt+tt- and D°D°Tr° /D°D*° modes. Three 
groups [13, [1^, [11] have devoted themselves to this study 
and the conclusions are quite different. One group [84] 
conclude that the X(3872) tends to be a virtual state of 
D^D*^, while another group's fit ^8§\ favors the loosely 
bound state explanation. Most recently, with the Belle's 
new data (zOjIzl], authors of Ref. [S^ gave a more system- 
atic study on this topic and found that there may need 
to be two near-threshold poles to account for the data, 
one from the D'^D*^ component and the other from the 
charmonium state x'ci ■ 

To sum up for the discussion in this subsection, we find 
that the Xci (2P)-dominated charmonium interpretation 
for the X(3872) may account for (i) the El transition 
rates to J/'0 and ijj{2S) and their ratio ((53)) : (ii) the large 
production rates in B decays and equal rates for 5+ and 
S"; (iii) the large production rate in p—p collisions at the 
Tevatron; (iv) the isospin violating decay to J/ijjp. More- 
over, in the screened potential model the mass of Xci(2P) 
is predicted to take a lower value than the quenched po- 
tential model. However, though the mass of Xci(2P) can 
be lowered by coupling to D^D*^ , one can not provide 
a quantitative explanation for the extreme closeness of 
X(3872) to the D^D*^ threshold (say, within 0.5 MeV), 
which is the most favorable motivation for the molecule 
interpretation. 



E. X(3940),X(4160) 

The X(3940) was found by the Belle Collaboration [1] 
in the recoiling spectrum of J/^ in the e"*"e~ annihi- 
lation process e^e~ — > J/if) + anything and e^e~ — > 
J/V-J + DD* . The later was studied further with higher 
statistics by Belle The mass and width of X(3940) 
are determined to be 

M(X(3940)) = 39421^ ±6 MeV, (55) 
r(X(3940)) = 37li^±8MeV. (56) 

Meanwhile, they also found the X(4160) in the D*D* 
mode in the process e+e" — > J /"ip + D*D* with a signif- 
icance of 5.1(7. The mass and width of the X(4160) are 
given by 

M(X(4160)) = 4156^20 ± 15 MeV, (57) 
r(X(4160)) = 139tg" ± 21 MeV. (58) 

Besides, there is a structure around 3880 MeV in the DD 
spectrum in e"'"e~ — > J/ip + DD. However, it is too wide 
to present a resonance shape sufficiently. 

Both of the two X-states have large production rates 
in these processes [s^j. This fact implies that the charge 
parities should be even since the charge odd state as- 
sociated with J/^ needs to be produced via two pho- 
ton fragmentation, which is expected to be highly sup- 
pressed [4^. On the other hand, the only known char- 
monium states that are produced in this way are 77c, 77^ 
and Xco fssj . and this double charmonium production 
phenomenon can be explained in the framework of non- 
relativistic QCD The production rates of X(3940) 
and X{A160) [13] are both as large as those of rjc, v'c 
XcO- This suggests that the two X-states could be either 
pseudoscalar like rjc or scalar like XcO (see Ref. [1^ for 
more detailed discussions). 

The observation that the dominant decay mode of 
X(3940) being DD* and the lack of evidence for the 
DD decay mode [1, [13] indicates that it can not be 
a scalar but can be a good candidate for the j]c{3S). 
The main problem is the low mass of X(3940) as the 
ric{3S). Although lower than that in the quenched po- 
tential model [5[ by 50 MeV or more, the mass of r]c{3S) 
in our screened potential model (see Table|T| is still larger 
than the observed mass (|55p by about 50 MeV. Moreover, 
the mass splitting between X(3940) and ^'(4040) is larger 
than that between rj'^ and -0', which looks quite unnatu- 
ral. But it may be due to the coupled channel effects [3], 
which will further lower the ric{3S) mass hopefully. 

The dominant decay mode of X(4160) is D*D* [13], 
and the other modes, such as DD and DD*, were not 
seen. Thus the charmonium candidates can be 77c (45) 
or Xcoi'iP), whose masses are 4250 MeV and 4131 MeV 
in our model prediction, respectively. Evidently, the 
mass of X(4160) in ([57]) favors the Xco(3-P) candidate. 
The Xco(3-P) can not decay into DD*, and the decay 
XcO (3^*) DD is expected to be strongly suppressed by 
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the form factor and the effects induced by the nodes of 
the 3P wave function, just hke the case of suppressed 
■0(4040) DD decay Q- The main problem of this as- 
signment may be why the Xco(2-P) state is not found in 
the similar process. One possible account [45] is that the 
broad peak around 3880 MeV in the DD spectrum jSTj 
mentioned above could be the missing Xco(2-P) state, 
since its mass in our model is about 3842 MeV and just 
lies within the bump (note, however, that this bump 
might not be a resonance |90l|l. In addition, the mea- 
surements on angular distributions can be used to test 
the two possible assignments, rjd^S) and Xco(3P), for 
the X(4160) [45i]. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we try to incorporate the color-screening 
(string breaking) effect due to light quark pair creation 
into the heavy quark-antiquark long-range confinement 
potential, and investigate the effects of screened po- 
tential on the spectrum of the charmonium especially 
the higher charmonium. We calculate the masses, elec- 
tromagnetic decays, and El transitions of charmonium 
states in the nonrelativistic screened potential model, 
and propose possible assignments for the newly discov- 
ered charmonium or charmonium- like states, i.e., the so- 
called "X, F, Z" mesons. We find that compared with 
the unscreened potential model, the masses predicted in 
the screened potential model are considerably lower for 
higher charmonium states. For example, the predicted 
Xc2(2P) mass well agrees with that of the Z(3930), and 
the mass of ij{5S) rather than ^p{4:S) is compatible with 



that of -0(4415). As a result of the compressed mass spec- 
trum in our model, most of the X, Y, Z states might be 
accomodated in the conventional higher charmonia. In 
particular, the discovered four Y states in the ISR pro- 
cess, i.e., r(4008),y(4260),r(4320/4360),r(4660) may 
be assigned as the -0(35'), 0(45*), (31?), 2p{QS) states re- 
spectively. The X(3940) and X(4160) found in the dou- 
ble charmonium production in e^e~ annihilation may 
be assigned as the rye (35*) and Xco(3P) states respec- 
tively. Based on the calculation of El transition widths 
for Xci(2P) 7 J/0 and Xci{2P) 70(25) and other 
results, we argue that the X(3872) may be a Xci(2-P) 
dominated charmonium state with some admixture of the 
£)0£)*o component. The problems encountered in these 
assignments and comparisons with other interpretations 
for these X, Y, Z mesons are discussed in detail. We em- 
phasize that more theoretical and experimental investi- 
gations are urgently needed to clarify these assignments 
and other interesting interpretations. In particular, we 
hope experiments at BESIII and SuperBelle in the fu- 
ture will be crucially useful in searching for new hadrons 
including charmonium-like states and testing the theo- 
retical interpretations. 
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TABLE I: Experimental and theoretical mass spectrum of charmonium states. The experimental masses are PDG averages. 
The masses are in units of MeV, while the averaged radii are in units of fm. The results of our screened potential model are 
shown in comparison with that of Ref.Q| including the NR and GI models[^. 





State 


Expt. 


Theor. 


of ours 


Theor. 


of Ref.[5] 








Mass 




NR 


GI 


IS J/iPil-'Si) 


3096.916 ±0.011 


3097 


0.41 


3090 


3098 


2980.3 ± 1.2 


2979 




2982 


2975 


2S 


V''(2^Si) 


3686.093 ± 0.034 


3673 


0.91 


3672 


3676 




3637 ± 4 


3623 




3630 


3623 


3S 


t/)f3^Si) 
rye (3"^ So) 


4039 ± 1 


4022 


1.38 


4072 


4100 






3991 




4043 


4064 


4S 


i/;f4^Si) 

r7.f4^Sn1 


4263+q 

— y 


4273 


1.87 


4406 


4450 






4250 




4384 


4425 


5S 


V'(5^Si) 


4421 ± 4 


4463 


2.39 










4446 








uo 


n ("fi-^Snl 




4608 


2.98 










4595 








1 p 




3556.20 ± 0.09 


3554 


0.71 


3556 


3550 






3510.66 ±0.07 


3510 




3505 


3510 






3414.75 ±0.31 


3433 




3424 


3445 




llc\-L ^ IJ 


3525.93 ±0.27 


3519 




3516 


3517 


2P 




3929 ± 5 ± 2 


3937 


1.19 


3972 


3979 




Yi f2^Pi ) 




3901 




3925 


3953 








3842 




3852 


3916 




/ic(2iPi) 




3908 




3934 


3956 




X2(3^P2) 




4208 


1.67 


4317 


4337 




Xi(3'Pi) 




4178 




4271 


4317 




Xo(33Po) 




4131 




4202 


4292 






4184 




4279 


4318 


ID 


V3(1'D3) 




3799 


0.96 


3806 


3849 




V2(l'D2) 




3798 




3800 


3838 




t/'(l'Di) 


3775.2 ±1.7 


3787 




3785 


3819 




77c2(l'D2) 




3796 




3799 


3837 


2D 


V'3(2-^D3) 




4103 


1.44 


4167 


4217 




V'2(23D2) 




4100 




4158 


4208 




V'(2='Di) 


4153 ± 3 


4089 




4142 


4194 




??c2(2^D2) 




4099 




4158 


4208 


3D 


V'3(3-*D3) 




4331 


1.94 








V2(3'D2) 

V(3='Di) 

77c2(3lD2) 




4327 
4317 












4326 









TABLE IL Leptonic widths (in units of KeV) for charmonium states without S-D mixing in the screened potential model. The 
widths calculated with and without QCD corrections are marked by Fee and respectively. The experimental values are 
taken from PDG 



state 


pU 

ee 


Tee 


^ ee 


l^S'i(3097) 


11.8 


6.60 


5.55 ±0.14 ±0.02 


2^51(3686) 


4.29 


2.40 


2.33 ±0.07 


3^ Si (4039) 


2.53 


1.42 


0.86 ±0.07 


4^ Si (4263) 


1.73 


0.97 




5^51(4421) 


1.25 


0.70 


0.58 ±0.07 


6-' Si (4664) 


0.88 


0.49 




l^i3i(3775) 


0.055 


0.031 


0.259 ±0.016 


2^i3i(4153) 


0.066 


0.037 


0.83 ±0.07 


3-'Di(4361) 


0.079 


0.044 
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TABLE III: Two-photon decay widths (in units of KeV) of pseudoscalar (^5*0), scalar (^Po), and tensor (^^2) charmonium states. 
Charmonium masses are in units of MeV. 







Theory 


PiYTi r>r 1 in PTil" 






Ref.[22l 


Ref.[33] 


Ref.f34] 


Ref.f35] 


Ref.f36] 


Ref.f37] 
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TABLE IV: El transition rates of charmonium states in the non-screened potential model [3| (marked by NR) and our 
screened potential model (those calculated by the zeroth-order wave functions are marked by SNRo and those by the first- 
order relativistically corrected wave functions are marked by SNRi). 
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